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Abstract. A non-relativistic Fokker-Planck equation for the electron distribution 
function is formulated incorporating the effects of stochastic acceleration by whistler- 
mode waves and Coulomb collisions. The stationary solution / to the equation, subject 
to a zero-flux boundary condition, is found to be a generalized Lorentzian (or kappa) 
distribution, which satisfies / oc t!~ 2<K+1 ) for large velocity v, where k is the spectral 
index. The parameter k depends strongly on the relative wave intensity R. Taking 
into account the critical energy required for resonance of electrons with whistlers, we 
calculate a range of values of R for each of a number of different space plasmas for 
which kappa distributions can be expected to be formed. This study is one of the first 
in the literature to provide a theoretical justification for the formation of generalized 
Lorentzian (or kappa) particle distribution functions in space plasmas. 
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1. Introduction 

In the natural space environment, e.g., planetary magnetospheres, and the solar 
wind (and actually many other kinds of astrophysical objects), plasmas are generally 
observed to possess a particle distribution function with a non-Maxwellian high-energy 
tail. The distribution function typically has a power-law tail in energy [distribution 
oc (particle energy) _ ( K+1 )] ) and frequently can be modeled by a generalized Lorentzian 
(kappa) distribution. For a variety of observational data, see Vasyliunas [1968], Gosling 
et al. [1981], Armstrong et al. [1983], Christon et al. [1988], Divine and Garrett 
[1983] and references contained in Summers and Thome [1991], Collier [1993], and 
Mace and Hellberg [1995]. Associated with the generalized Lorentzian distribution has 
been developed the modified plasma dispersion function which is an effective tool for 
investigating waves and microinstabilities in space plasmas \\Summers and Thome, 
1991| ; \Mace and Hellberg, 1995| ; \Summers et al., 1994| ; |1996|| . Physical mechanisms that 
can produce power-law particle distributions include stochastic acceleration by plasma 
wave turbulence, and collisionless shocks. Kappa-like distribution functions which are 
Maxwellian at low energies and power law at high energies can also be produced by 
velocity-space Levy flight probability distributions [e.g., Collier, 1993]. It is the process 
of stochastic acceleration by plasma wave turbulence with which we shall be concerned 
in this paper. Many researchers have studied stochastic acceleration of charged particles 
by various modes of plasma wave, e.g., Gurevich [1960], Kennel and Engelmann [1966], 
Melrose [1980, 1986], Steinacker and Miller [1992], and Schlickeiser [1997]. Hasegawa 
et al. [1985] showed that superthermal radiation can enhance velocity-space diffusion 
so as to produce a power-law distribution. However, these authors did not estimate 
values of the intensity of the radiation field and other plasma parameters necessary to 
obtain a reasonable value for the power law. The source of the free energy to excite 
such intense radiation was another question left unaddressed. Melrose [1980] pointed 
out that the equilibrium distribution is a power law only if the momentum diffusion 
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coefficient D(v) is proportional to v _1 , where v is the particle velocity. For example, 
for Langmuir waves, D(v) oc i> _1 ; but this requires the turbulent spectrum I{k) to be 
independent of the wave number k, which is an unlikely possibility [Melrose, 1980]. 
Stochastic acceleration by plasma waves is characteristically slow due to the diffusive 
nature of the scattering process. Nevertheless, the process may be applicable, for 
instance, in planetary magnetospheres in cases when typical electron acceleration times 
are of the order of tens of hours. We note that Schlickeiser [1997] has recently obtained 
a diffusion coefficient varying as D(v) oc v _1 when the charged particles are accelerated 
by whistler-mode waves. The momentum diffusion coefficient results from the resonant 
interaction between particles and whistler-electron cyclotron plasma waves. Similar 
formulae have also been obtained by Dermer et al. [1996]. Numerical calculations 
of the momentum diffusion coefficients for particle-whistler interaction were given by 
Steinacker and Miller [1992] and Pryadka and Petrosian [1997], but no analytical 
formulae were obtained for medium energetic electrons (10 keV to a few hundreds 
of keV), which limits their application by other authors. There are several physical 
processes that can drive the whistler instability. For example, anisotropies in the plasma 
velocity distribution can make whistlers unstable. Moreover, a kappa distribution 
can itself enhance the whistler-mode instability, as compared with the Maxwellian 
distribution [e.g., see Summers and Thome, 1992; Xue et al, 1993; Mace, 1998]. 

In the following section, we shall develop a non-relativistic theory for the formation 
of kappa distributions by means of whistler-mode wave stochastic acceleration. 

2. Theory 

We assume a Kolmogorov-type magnetic turbulence power spectrum in the wave 
number for the whistler-mode waves, namely, 

I = W 9 > (!) 



5 



with q > 1, where k\\ > k m i n , k\\ is the parallel wave number, and Iq = (q — l)(5i?) 2 fc^ is 
the energy density of the whistler turbulence. It is natural to assume that k min VA = 
for whistler waves, where Va and Qi are respectively the Alfven speed and the 
non-relativistic ion gyrofrequency. In this paper we are concerned with the situation 
in which pitch-angle anisotropy is small and the particle distribution is isotropic. For 
higher energy electrons, it is generally the case that the time-scale for pitch angle 
diffusion is much less than that for momentum diffusion, and so the assumption of 
isotropy is justified. For lower energy electrons, where the pitch angle and momentum 
diffusion rates are comparable, we assume that isotropization by Coulomb collisions 
occurs. The momentum diffusion coefficient can then be averaged over the pitch angle. 
Herein, we employ the diffusion coefficient given by Schlickeiser [1997], namely, 



ir(q - 1) [5B\ 2 2+q ( k min c x ' ! i 



where Q e is the non-relativistic electron gyrofrequency, Pa = Va/c, /3 — v/c, c is the 
speed of light, and J\y is a weakly varying function of v. Henceforth in this paper, the 
Kolmogorov turbulent spectrum (q = 5/3) is adopted. It then follows that Jw is of 
order unity. In accordance with the whistler wave dispersion relation and the resonance 
condition \\Melrose, 1986|| , we find that electrons can resonate with whistlers provided 
that (3 > (m p /m e ) l l 2 (3a- The diffusion coefficient (0) is valid in the region where 
(m p /m e ) 1//2 Pa < P < (fn P /me)PA- However, lower-energy electrons can be accelerated 
by high wave-number whistlers as the electron cyclotron branches of the dispersion 
relation are approached \\Petrosian, 1994 . 

Coulomb interactions result in both momentum diffusion and friction in the kinetic 
equation for the particle distribution function. As we are assuming the distribution to 



be isotropic, we can ignore the angular contribution to Coulomb collisions \Hinton, 



1983], giving the equation 

dt v 2 dv \ CJ ' 2 " c dv 



v - 1» { FJ + i D , c m , (3) 



6 



where the frictional coefficient F c is given by 

, D\\ c -D± c , ldD\\ c 
v 2 dv 

v a is the slowing-down rate, and D^ c and Dn c are the respective perpendicular and 
parallel diffusion coefficients. When the velocity of the particle v is much larger than 
the thermal velocities of the electrons and ions [v e <C v,Vi <C v,v Bi i = (2T e> i/m e) i) 1 ' 2 ], we 
have by Hinton [1983], 

v s = 3n e -|, (5) 
tr 

D ±c = (6) 



D = n e T e (vl + vf) _ n e T e v 2 e 



where n e is the electron number density, 



4vre 4 log 10 A 
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and log 10 A ~ 20 is the Coulomb logarithm. Finally, modifying Equation (JJ) to include 
acceleration by whistler-mode waves, we obtain the following Fokker-Planck equation to 
describe the evolution of the distribution function: 

df _ 1 d v 2f F j + 1 D d f\ + 1 9 ( v z D d f\ (g) 
dt v 2 dv \ c 2 " c dv J v 2 dv I dv J ' 

where Z) is given by Equation (0). With a zero- flux boundary condition, the stationary 
solution of Equation (^) is given by 

*-*->{- f^h**}- (9) 

Substituting the expression (2) for D in Equation @, we obtain 
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-(k+1) 



.4 



1 + 



K0 2 



(10) 



where 



K+ 1 



24e 4 n e log 10 A 



11^ 



Q e m p m e c 3 Ri3\Jw ' 

A is a constant of integration, i? = (5B/B ) 2 is the power of the wave turbulence, and 
= [(/c + l)/^] 1 / 2 ^. Normalizing the distribution using Airv 2 f dv = 1 gives 

,_ 1 r(« + i) 

(itk) 3 / 2 6 3 T(k- 1/2) 

where T is the gamma function. The distribution ([H]) is formally a generalized 
Lorentzian (kappa) distribution, with the spectral index k given by Equation ([TTD . 
Parenthetically, it can be noted that setting D = in Equation ([]) leads to the recovery 
of the Maxwellian distribution / = Aexp(— v 2 /v 2 ). From Equation ( |IT]) we determine 
that k is given by 

« + 1 = 0.9 x 10~ 25 n 5 J 2 Bq 4 R- 1 J w\ (12) 
where n e is in cm -3 , and B in Gauss. 



3. Discussion 

In the foregoing, we have determined that stochastic acceleration by whistler-mode 
turbulence can produce a particle distribution with a power-law tail, i.e., / oc t>~ 2 ( K+1 ). 
By Equation fll2|) , we note that the value of the spectral index k depends strongly on 
the relative wave intensity R. In Table [l], for various space plasmas, we have listed 
typical values of the parameters n e ,B ,PA and the thermal energy(temperature). The 
basic plasma parameters in this table for the solar wind were taken from Kivelson 
and Russell [1995], for the Earth's plasma sheet from Christon et al. [1988], for the 
Jovian magnetosphere from Dessler [1983] and Divine and Garrett [1983], and for 
the Saturnian magnetosphere from Sittler et al. [1983]. Because only the energetic 
electrons with velocities satisfying (3 > (m p /m e ) l l 2 (3a can resonate with first-harmonic 



Table HI 



8 



whistlers, we list in Table 1 the critical energy E c corresponding to the critical velocity 

P c = (nip/me) 1 / 2 Pa- If the critical energy E c is much larger than the electron thermal 

energy, then the whistlers can only accelerate non-thermal high energy particles, 

for which the Coulomb momentum diffusion coefficient D» c is much less than the 

acceleration diffusion coefficient. In this case, an injection mechanism is required to 

pre-accelerate the particle to the required critical energy E c . It is straightforward 

to show that this situation applies, for instance, to solar flares. In fact, there is no 

observational evidence that electrons have a kappa distribution in solar flares. If 

the critical energy E c is of the same order as the thermal energy, then whistlers can 

accelerate the thermal particles into the high-energy tail of the distribution, and so a 

kappa distribution can be produced. This is essentially the case for all the space plasmas 

shown in Figure |I|. For a specific plasma characterized by its electron number density Figure [I] 

n e and background magnetic field Bo, formula (|12|) gives the spectral index k produced 

by whistler turbulence of power R. In Figure |], for each of the given space plasmas, we 

show the range of .R-values associated with the range of K-values that might be typically 

observed, namely 2 < k < 8. In addition, we also show the value of the critical energy 

E c for each plasma. The results imply that relatively weak whistler-mode turbulence 

could produce power-law spectra in the Earth's plasma sheet and in the Jovian and 

Saturnian magnetospheres, whereas stronger turbulence is required in the solar wind. 

In conclusion, we should point out both the merits and the limitations of this 
study. Of course, it is not the case that all electron power-law spectra are produced 
by whistler-mode turbulence, or that whistler-mode turbulence will always give rise to 
electron power-law spectra. Nevertheless, under certain restrictions, we have shown that 
the production of electron power-law energy spectra by stochastic acceleration due to 
whistlers is a viable possibility. Apart from the non-relativistic fromework, limitations 
of the study include the assumptions of isotropy, a zero-flux boundary condition, and 
a momentum diffusion coefficient based on wave propagation strictly parallel to the 
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background magnetic field. The intrinsic value of the study is that it is one of the first 
in the literature to present a firm theoretical basis for the formation of generalized 
Lorentzian (or kappa) particle distribution functions in space plasmas. 
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Figure 1. The range of values of the relative wave intensity R , calculated from Equation 
((Hp, corresponding to the range of k- values, 2 < k < 8, for space plasmas for which 
typical parameters and E c - values are given in Table 1. 
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Table 1. Typical parameters for different space plasmas. Also given are values for the 
critical energy E c which is the minimum particle energy required for particle acceleration by 
whistler-mode waves. 

Plasma parameters Solar wind a Earth's plasma sheet b Jupiter c Saturn d 



Density n e (cm 3 ) 7 

Magnetic field B (G) 7 10" 5 

Alfven speed parameter (3a 1.9 10" 
Thermal energy (eV) 12 



1 2 10 3 0.14 

2 10" 4 1.8 1(T 2 6 1(T 5 

1.5 lCT 3 2.8 10- 3 1.1 1CT 3 
1000 1000 300 



Critical energy E c (eV) 17 1500 3500 500 



a At 1 AU, A.J. Hundhausen, in Kivelson and Russell [1995], page 92. 
b At L=14, Christon et al. [1988]. 

c At L=6, e.g., Dessler [1983] and Divine and Garrett [1983]. 
d At L=15, Sittler et al. [1983]. 
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